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Abstract 

The genetic structure of populations of closely related, sympatric species may 
hold the signature of the geographical mode of the speciation process. In fully 
allopatric speciation, it is expected that genetic differentiation between species 
is homogeneously distributed across the genome. In nonallopatric speciation, 
the genomes may remain undifferentiated to a large extent. In this article, we 
analyzed the genetic structure of five sympatric species from the plant genus 
Heliotropium in the Atacama Desert. We used amplified fragment length poly- 
morphisms (AFLPs) to characterize the genetic structure of these species and 
evaluate their genetic differentiation as well as the number of loci subject to 
positive selection using divergence outlier analysis (DOA). The five species form 
distinguishable groups in the genetic space, with zones of overlap, indicating 
that they are possibly not completely isolated. Among-species differentiation 
accounts for 35% of the total genetic differentiation (F ST = 0.35), and F ST 
between species pairs is positively correlated with phylogenetic distance. DOA 
suggests that few loci are subject to positive selection, which is in line with a 
scenario of nonallopatric speciation. These results support the idea that sympat- 
ric species of Heliotropium sect. Cochranea are under an ongoing speciation 
process, characterized by a fluctuation of population ranges in response to 
pulses of arid and humid periods during Quaternary times. 
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Introduction 

Based on the spatial segregation of the involved popula- 
tions, three major geographical categories have tradition- 
ally been recognized for the process of speciation. In the 
case of allopatric speciation, geographically isolated popu- 
lations differentiate into distinct species and the specia- 
tion event occurs in the absence of gene flow. Sympatric 
speciation, on the other hand, is a process that occurs in 
a single (geographical) population and in the presence of 
gene flow. Besides these two extreme forms of spatial seg- 
regation, partial segregation of the diverging populations 
can occur. The diversification of such populations with 
adjacent geographical ranges into distinct species, despite 



limited interbreeding at their contact zone, characterizes 
the process of parapatric speciation. 

Classification of the speciation process into the dis- 
crete allopatric, parapatric, and sympatric modes is now 
considered an oversimplification by several authors who 
have argued that speciation studies should focus on mod- 
eling and estimating parameters that describe the diversi- 
fication process rather than on classifying the speciation 
mode according to geographical criteria (e.g., Wu 2001; 
Butlin et al. 2008; Fitzpatrick et al. 2008; Coyne 2011). 
Nevertheless, theoretical models have illustrated the possi- 
bility of sympatric speciation (Bolnick and Fitzpatrick 
2007), and its occurrence in nature has been shown in 
several studies (e.g., Barluenga et al. 2006; Savolainen 
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et al. 2006; Papadopulos et al. 2011). The relative scarcity 
of reports on the occurrence of speciation with gene flow 
in nature is thereby likely due to the fact that allopatric 
speciation is commonly accepted as the null model 
(Coyne and Orr 2004) and nonallopatric alternatives are 
usually not considered in speciation studies (Johannesson 
2010). Despite this apparent bias in the treatment of 
allopatric and sympatric speciation models, the results of 
Papadopulos et al. (2011) show that sympatric speciation 
may be more common than previously accepted and 
highlight the relevance of gene flow between incipient 
species (Wu 2001; Fitzpatrick et al. 2009; Coyne 2011). In 
order to comply with the current views on the different 
modes of speciation, in this study, we refer to sympatric 
speciation as a speciation process that started in the pres- 
ence of panmixia. 

Sympatry or co-occurrence of closely related species can 
either result from a sympatric speciation process or from 
secondary contact due to range expansion after speciation. 
Theoretical models as well as empirical studies have shown 
that the genetic structure of closely related sympatric 
species holds a signature of the geographical mode of the 
speciation process (e.g., Charlesworth et al. 1997; Wilding 
et al. 2001; Savolainen et al. 2006; Nosil et al. 2009; Stras- 
burg et al. 2012; Via 2012). Under the allopatric scenario, 
genetic variation tends to be uniform across the genome 
due to a large proportion of the genome changing through 
a combination of divergent selection, differential response 
to similar selective pressures and genetic drift (Nosil et al. 
2009; Strasburg et al. 2012). In contrast, in the extreme 
case of sympatric speciation, gene flow between the incipi- 
ent species can homogenize most of the genome, except 
for the loci that experience strong divergent selection pres- 
sures or regions that are tightly linked with these loci 
(Charlesworth et al. 1997; Strasburg et al. 2012; Via 2012; 
Via et al. 2012). Closely related species currently living in 
local sympatry thus provide an opportunity to examine 
these models empirically and to infer possible modes of 
speciation. However, closely related species tend to have 
similar ecological characteristics, and their sympatry is 
constrained by competition (HilleRisLambers et al. 2012). 
For this reason, cases of clades containing multiple sym- 
patric and closely related species are rare. 

In the Atacama Desert, there are at least two docu- 
mented groups of closely related plant species that are to 
a large extent sympatric: Nolana (Dillon et al. 2009) and 
Heliotropium (Luebert and Wen 2008; Luebert 2013). All 
but one of the Heliotropium species that inhabit the Ata- 
cama Desert belong to the monophyletic section Cochra- 
nea (Miers) I.M.Johnst. (Luebert and Wen 2008). This is 
one of the most speciose groups of the Atacama Desert, 
with all of its 17 species inhabiting that area (Luebert 
2013). Phylogenetic studies are in line with the idea that 



Heliotropium sect. Cochranea underwent an adaptive radi- 
ation toward the early Pliocene in response to the aridifi- 
cation of the Atacama Desert (Luebert and Wen 2008; 
Luebert et al. 2011a,b). However, sister relationships are 
largely unresolved, and little is known about the ecologi- 
cal and genetic forces as well as the spatial and temporal 
dynamics implied in its diversification. Incipient specia- 
tion is likely common in this group. Based on the ideas 
of Stebbins (1952) for speciation in arid zones, Luebert 
and Wen (2008) postulated that late Neogene and Qua- 
ternary diversification of Heliotropium sect. Cochranea 
may have occurred along with temporal pulses of spatial 
contraction and expansion of populations, generating 
alternating periods of isolation and reunion of diverging 
populations. This model of diversification implies allopat- 
ric speciation, because genetic isolation would occur 
during the contraction phases, and is consistent with 
long-term fluctuations of aridity in the Atacama Desert 
(Latorre et al. 2006). In fact, most evolutionary studies 
on Atacama Desert plant groups suggest a predominance 
of allopatric speciation (Gengler-Nowak 2002; Luebert 
et al. 2009; Viruel et al. 2012). Present distributions in 
Heliotropium sect. Cochranea show two major centers of 
diversity with, respectively, eight and six of 17 species 
sharing the same geographical area, and up to 6 species 
can occur in local sympatry within these areas. If the 
model postulated by Luebert and Wen (2008) holds, 
sympatry would have been achieved through secondary 
contact. Under this scenario, we would expect absence of 
gene flow among sympatric species associated with uni- 
form genetic divergence across the genome. Under a sce- 
nario of sympatric speciation, however, we expect to find 
genetically differentiated populations with zones of 
genetic contact between species reflecting various levels of 
gene flow. Genetic divergence between species is then 
expected to occur in only a few loci that are subject to 
selection. However, depending on the time frame of 
allopatry, the pattern of gene flow may be very difficult 
to distinguish from the pattern expected under a scenario 
of speciation with gene flow (Strasburg et al. 2012). 

In this study, we evaluate the genetic differentiation 
among five sympatric and closely related species of the 
plant genus Heliotropium L. (Heliotropiaceae) in the 
Atacama Desert using amplified fragment length polymor- 
phism (AFLP) data. The following questions are addressed: 

1 Are closely related and sympatric species of Heliotropi- 
um genetically isolated? 

2 Is phylogenetic relatedness of Heliotropium species 
associated with gene flow? 

3 Can a signature of allopatric, sympatric, or ongoing 
divergence with secondary gene flow in sympatry be 
found in the genetic structure of the Heliotropium 
species? 
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Material and Methods 

Plant material and DNA extraction 

The main study site is located in the area around the 
village of Totoral (S27°54'5", W70°57'34") in the coastal 
Atacama Desert of northern Chile. For population 
genetic analyses, we aimed at collecting leaves of at least 
15 specimens and one voucher herbarium specimen of 
each of the following Heliotropium species: Heliotropium 
filifolium (Miers) I.M.Johnst., H. floridum (A.DC.) Clos, 
H. longistylum Phil., H. megalanthum I.M.Johnst., and 
H. sinuatum (Miers) I.M.Johnst. These are morphologi- 
cally defined entities according to the species concepts of 
Luebert (2013). A total of five individuals of these spe- 
cies from the neighbor locality of Carrizal Bajo (S28°4' 
54", W71°8'49") was also included in the study. Samples 
were collected in the austral spring in 2003, 2005, and 
2011 (see Table SI), and leaf material was dried and 
stored in silica gel until further processing. DNA was 
isolated using the Nucleospin Plant® II Kit (Macherey- 
Nagel GmbH, Diiren, Germany) following manufacturer 
instructions. 

Amplified fragment length polymorphism 

Amplified fragment length polymorphism was used to 
evaluate genetic differentiation within and between the 
species under study. AFLP analysis followed the protocol 
of Vos et al. (1995). Restriction enzymes EcoRl and Msel 
were employed for digesting extracted DNA, and pream- 
plification was performed using oligonucleotide primers 
£coRI + A / Msel + C. Three different primer combina- 
tions were used for the selective amplification: EcoRl + 
ATG / Msel + CAG, £coRI + ATA / Msel + CAC, and 
£coRI + ATT / Msel + CAA. The selective EcoRl primers 
were labeled at their 5'-end with 6-carboxyfluorescein 
(6-FAM), and the selective PCR products were separated 
by electrophoresis using the Applied Biosystems DNA 
Analyzer 3730 (fragment analysis was performed by GATC 
Biotech AG (Konstanz, Germany)). Fragment length esti- 
mation and allele scoring were performed in GeneMarker® 
vl.95 (Softgenetics, State College, PA) using the GS500 
size standard and the default settings of the AFLP analysis 
parameters, except for the minimum intensity of the peak 
detection threshold, which was increased from 100 to 500. 
AFLP fragments of the same size generated in different 
specimens were assumed to be homologous, and the rela- 
tive intensities of the fragments were not considered. Loci 
for which the minor allele frequency (MAF) was at least 
3% were considered polymorphic, and only the polymor- 
phic loci were retained for further analysis. Reproducibil- 
ity of the AFLP method was assessed by comparing 



replicate analyses of three specimens of each of the species 
H. longistylum, H. filifolium, and H. megalanthum. 

Data analysis 

The level of linkage disequilibrium (LD) between AFLP 
markers and its significance was assessed by calculating r 2 
with the package LDcorSV vl.2 in R v2.15.2 (R Develop- 
ment Core Team 2013). An exact test of LD was 
performed in Genepop v4.2 (No. of steps in Markov 
chain = 10,000, No. of dememorization steps = 1000; 
Rousset 2008), using the subset of polymorphic AFLP 
markers in each species separately and a Bonferroni 
correction to obtain a global significance of 0.05. Genetic 
diversity within species was estimated as the percentage of 
polymorphic AFLP markers and was corrected for unequal 
sample sizes using the method of multiple random sub- 
sampling with a sample size of seven and 50 subsamples 
(Leberg 2002). Genetic distances among individuals were 
calculated as pairwise asymmetric binary distances or Jac- 
card dissimilarities between individual samples (see e.g., 
Kosman and Leonard 2005) and visualized by multidi- 
mensional scaling (MDS) and hierarchical clustering 
(complete linkage) in R v2.15.2 with packages adegenet 
vl.3-8 (Jombart 2008) and ade4 vl.5-2 (Dray et al. 2007). 

In order to evaluate the association between phyloge- 
netic relatedness and gene flow, pairwise phylogenetic 
distances between species were calculated with the R pack- 
age APE v3.0-8 (Paradis et al. 2004). A phylogenetic tree 
was generated by reanalyzing the sequence data obtained 
by Luebert and Wen (2008). These data consist of DNA 
sequences of three plastid regions (ndhE, rpsl6, trnL-trnF) 
and the nuclear ribosomal ITS region. In this dataset, each 
species of Heliotropium sect. Cochranea, except for 
H. megalanthum, H. myosotifolium (A.DC) Reiche, and 
H. krauseanum Fedde, is represented by two or more sam- 
ples from different localities. For one additional individual 
of each of these three species, additional data (GenBank 
accessions Nrs. KF301622-KF301628) were generated by 
sequencing the same four markers following the protocols 
described by Luebert and Wen (2008). Information of 
sequences used in this study for phylogenetic analysis is 
detailed in Table S2. Sequence data were analyzed using a 
species-tree approach implemented in the software 
* BEAST v. 1.6.1 (Heled and Drummond 2010). Each locus 
was analyzed using an uncorrelated lognormal relaxed 
clock model. Substitution models were specified according 
to results obtained from analyses with Modeltest v.3.7 
(Posada and Crandall 1998). Modeltest uses the Akaike 
information criterion to select one among 24 possible sub- 
stitution models with different numbers of parameters. 
Likelihoods of the alignment given a tree are calculated 
under each possible model using neighbor-joining trees. 
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The final species tree was deposited in TreeBase (accession 
14422). 

Genetic differentiation within and among species was 
evaluated by analysis of molecular variance (AMOVA; 
Excoffier et al. 1992) using the R package pegas vO.4-4 
(Paradis 2010) with pairwise Euclidean distances between 
the individual AFLP phenotypes and a permutation proce- 
dure (1000 permutations) to assess the significance of the 
different variance components. In our model, AMOVA 
partitioned the total genotypic variance into components 
due to differences between species and differences between 
specimens within species. Gene flow between species was 
estimated using the infinite island model as the effective 
number of migrants per generation N e m = (1— _F ST )/(4P ST ) 
(Wright 1949). 

In order to detect loci under selection, divergence out- 
lier analysis (DOA) was conducted with the software 
Mcheza (Antao and Beaumont 2011). Pairwise compari- 
son of P ST values between the species involved in this 
study was conducted with default settings. The subsample 
size and the expected number of total populations were, 
respectively, set to 0 and 2 for each species combination. 
The expected Fst distribution under the null hypothesis 
of neutrality was calculated with 50,000 simulations and a 
false discovery rate of 0.1 within the 95% confidence 
interval in 10 different simulations. As calculation of the 
initial mean F S t found in the dataset includes all poten- 
tially selected loci (i.e., loci falling outside of the confi- 
dence interval), those loci were removed from calculation 
of the "neutral" mean F ST . The mean simulated f ST was 
forced to approximate the initial mean F sx by application 
of a bisection algorithm over repeated simulations. 

Results 

Amplified fragment length polymorphism analysis of 96 
specimens belonging to five different Heliotropium species 
with three different selective PCR amplifications generated 
a total of 311 scorable loci, 287 (92%) of which were 
polymorphic across the whole set of samples. Per selective 
PCR, the number of polymorphic loci varied between 82 
and 118 (with an average of 96; see Table 1). Reproduc- 
ibility of the AFLP method was high, with an average 
allele call correspondence between replicates of 97%. 
Based on the set of 287 polymorphic loci, 96 different 
genotypes could be distinguished and the number of 
scored AFLP fragments per specimen varied between 72 
and 154 (average of 120). Overall, linkage disequilibrium 
between AFLP markers was found to be low, with the 
average r 2 value within species varying between 0.054 
(H. megalanthum) and 0.17 (H. floridum; see Table 2) 
and the fraction of the marker pairs in significant linkage 
disequilibrium not exceeding 0.4% (data not shown). The 



Table 1. Overview of oligonucleotide primer combinations used for 
the selective PCR amplification in the AFLP analysis of 96 specimens 
belonging to five different Heliotropium species and the numbers of 
polymorphic loci that were generated for each primer combination. 



fcoRI + 


Mse\ + 


# Polymorphic 
markers 


ATG 


CAG 


87 


ATA 


CAC 


82 


ATT 


CAA 


118 


Total 




287 



percentage of polymorphic markers per species, corrected 
for unequal sample size, varied between 55% (H. longisty- 
lum, H. filifolium, and H. megalanthum) and 68% 
(H. floridum; see Table 2). 

Considering both the results of the hierarchical cluster 
analysis (Fig. 1A) and of the MDS (Fig. IB), all species 
included in this study tended to form well-defined genetic 
groups. In some cases, however, specimens did not cluster 
according to their presumed taxonomic identity: one 
H. filifolium specimen clustered together with H. sinua- 
tum and one H. floridum specimen clustered together 
with H. megalanthum. 

Bayesian phylogenetic inference was conducted using 
the substitution models selected from Modeltest for each 
molecular marker: ITS (GTR+I+T), trnL-trnF (HKY+T), 
rpsl6 (GTR+D, ndhF (GTR+I+r). The resulting species 
tree (Fig. ID) suggests that Heliotropium pycnophyllum 
and H. krauseanum are successive sister species to the 
remaining species of section Cochranea. Among the 
remaining species, H. filifolium and H. glutinosum appear 
in an unsupported sister relationship, while the other spe- 
cies form a well-supported clade, but with unresolved 
internal relationships. Most sympatric species in the study 
area fall in the latter, unresolved clade. 

Estimated migration rates declined with increasing phy- 
logenetic distance. Phylogenetic distances and N e m values 
among species showed a negative correlation (p = —0.87, 
P = 0.0009; see Fig. 1C and Table S3). The same pattern 
was observed when comparing the number of private loci 
between species (data not shown). 

Analysis of molecular variance indicated that approxi- 
mately two-thirds (65%) of the total variability of the 
data can be explained by genetic variation among speci- 
mens within the investigated species, while 35% of the 
total variability is due to genetic differentiation between 
species (Table 3). 

Divergence outlier analysis with Mcheza showed that a 
variable but limited number of loci putatively under posi- 
tive selection are detected in all pairwise comparisons 
of species (Table 4, Fig. SI). The distributions of frequen- 
cies of single-locus f ST revealed a similar pattern. When 
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Table 2. Genetic diversities, estimated as percentages of polymorphic AFLP markers and corrected for unequal sample sizes by multiple random 
subsampling (subsample size = 7; Leberg 2002), in five Heliotropium species and average pairwise linkage disequilibria between polymorphic 
markers within each species measured as the square of the correlation coefficient between allele frequencies (r 2 ); N, sample size; SD, standard 
deviation. 



% Polymorphic markers 



Species 


N 


Uncorrected 


Corrected (SD) 


^(SD) 


P 


Heliotropium longistylum 


17 


71 


55 (3) 


0.073 (0.12) 


0.0000 


Heliotropium filifolium 


25 


77 


55 (5) 


0.065 (0.14) 


0.0042 


Heliotropium megalanthum 


21 


72 


55 (3) 


0.054 (0.085) 


0.0000 


Heliotropium sinuatum 


26 


85 


60 (7) 


0.063 (0.11) 


0.0000 


Heliotropium floridum 


7 


68 


NA 


0.17 (0.21) 


0.0000 



(A) 



Height 

9 11 



(D) 



(B) 



CD O 
CO o 




(C) 



* . • 

• 4 

0 


• 


>>.*' 


e°° 


o ° 










1 

-0.1 


I I 
0.0 0.1 0.2 

PCo1(19.8%) 


0.3 


0 

P 


= 0.87, P= 0.0009 






8 


§ 


I 

0.005 


I I 
0.010 0.015 





0.98 



0.97 



0.99 



Phylogenetic distance 




H. chenopodiaceum 
H. longistylum 
H. sinuatum 
H. linariifolium 
H. philippianum 
H. taltalense 
H. inconspicuum 
H. stenophyllum 
H. floridum 
H. megalanthum 
H. eremogenum 
H. myosotifolium 
H. glutinosum 
H. filifolium 
H. krauseanum 

H. pycnophyllum 
OUTGROUPS 



H. filifolium was compared with H. floridum, H. longisty- 
lum, and H. megalanthum, a more even distribution could 
be observed (Fig. 2), indicating a higher level of differenti- 



Figure 1. (A) Dendrogram resulted from 
cluster analysis. (B) Scatter plot of the two first 
axis (percentage of variance explained 
indicated in brackets) of the multidimensional 
scaling. (C) Relation between interspecific 
phylogenetic distance and N e m between 
populations, depicting the curve adjusted with 
a linear model and indicating the Spearman 
correlation and its P-value. (D) Phylogenetic 
analysis with Bayesian Posterior probabilities 
>50 above branches, Heliotropium sect. 
Cochranea highlighted in gray, and studied 
species in bold and color-coded. Color codes 
for species are the same across figures. Black 
arrows (A-B) indicate individuals of one species 
intermixed with individuals of other species. 

ation (and therefore less gene flow). This is to be expected 
as H. filifolium is phylogenetically less related to the other 
species than the other species are to each other (Fig. ID). 
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Table 3. Results of an AMOVA indicating sums of square deviations (SSD), mean square deviations (MSD), number of degrees of freedom, and 
variance components. 





SSD 


MSD 


df 


Variance components 

a 1 % variation 


P-value 


Populations 


1539.909 


384.97737 


4 


18.862 


35.37 


0 


Error 


3135.882 


34.46024 


91 


34.460 


64.63 




Total 


4675.792 


49.21886 


95 


53.322 







Discussion 

Our results reveal a number of previously unknown find- 
ings regarding the five sympatric species of Heliotropium 
sect. Cochranea included in this study. First, the obtained 
phylogenetic relationships are well supported and consis- 
tent with previous works (Luebert and Wen 2008; Luebert 
et al. 2011a,b). However, the present analysis resolves a 
well-supported sister relationship between H. krauseanum 
and all other species except H. pycnophyllum (which is 
sister to all other species of this group). This relationship 
was not resolved in previous analyses (Luebert and Wen 
2008; Luebert et al. 2011a). Second, the five sympatric 
species seem to be genetically well differentiated from 
each other in spite of their relatively recent origin (Lue- 
bert and Wen 2008; Luebert et al. 2011b): They occupy 
distinct regions in the genetic space and, with the excep- 
tion of H. floridum versus H. megalanthum (F S t = 
0.1004), all species-pair comparisons have F ST estimates 
larger than 0.17 (see Table 4). This is reinforced by the 
fact that a considerable proportion of the total genetic 
differentiation (-35%) is due to among-species variation. 
Figure IB shows that genetic distance tends to be highest 
between H. filifolium and the other Heliotropium species, 
which is consistent with its phylogenetic position outside 
of a poorly resolved clade in which the remaining species 
are located (Fig. ID). However, sister relationships within 
this clade remain unknown. 

A few individuals that were morphologically assigned to 
a certain species clustered together with individuals of 



another species. Barring identification errors, this may 
indicate that, despite genetic differentiation among species, 
gene flow still occurs between the involved species. In fact, 
measures of gene flow calculated according to the finite 
island model also suggest that these species experience dif- 
ferent levels of gene flow and that levels of gene flow are 
correlated with phylogenetic relatedness. This picture is 
consistent with personal field observations on these sym- 
patric species that indicate clear morphological differentia- 
tion (Luebert 2013). Differentiation might be associated 
with different strategies to tolerate drought (Luebert et al. 
2011a) and how species are pollinated (Luebert 2013). In 
the latter sense, relative position and size of reproductive 
floral structures may be associated with different pollina- 
tors, a hypothesis that remains to be tested. Additionally, 
hybrids may occur in the areas of sympatry. In support of 
this, morphologically intermediate individuals between 
H. sinuatum and H. longistylum, between H. longistylum 
and H. floridum, and between H. floridum and H. maega- 
lanthum have been observed (Luebert 2013; pers. obs.). 
Field annotations are also consistent with the observed 
patterns of genetic differentiation. For example, the sam- 
ple Fili7 was annotated in the field to be a putative hybrid 
between H. filifolium and H. sinuatum, because it exhib- 
ited floral morphological characters of the former species, 
but leaf characters transitional to the latter. Our AFLP 
data suggest that this individual falls within the genetic 
space of H. sinuatum (see Fig. 1A). 

Comparisons of pairs of sympatric species with DOA 
indicate that few of the investigated loci are expected to 



Table 4. Results of the divergence outlier analysis for each species pair. Percentage of outliers as calculated in relation to all polymorphic loci. 



Species pair 


Dataset F ST 


"Neutral" F ST 


W outliers 


% outliers 


Heliotropium filifolium-Heliotropium floridum 


0.2489 


0.205 ± 0.018 


6 (0-9) 


2.28 


Heliotropium filifolium-Heliotropium megalanthum 


0.3184 


0.284 ± 0.016 


8 (0-14) 


3.04 


Heliotropium filifolium-Heliotropium longistylum 


0.3161 


0.281 ± 0.034 


4 (3-16) 


1.52 


Heliotropium filifolium-Heliotropium sinuatum 


0.2791 


0.256 ± 0.008 


7 (7-8) 


2.66 


Heliotropium floridum-Heliotropium megalanthum 


0.1004 


0.086 ± 0.008 


5 (5-7) 


1.90 


Heliotropium floridum-Heliotropium longistylum 


0.1771 


0.151 ± 0.009 


4 (4-7) 


1.52 


Heliotropium floridum-Heliotropium sinuatum 


0.1729 


0.137 ± 0.008 


8 (8-12) 


3.04 


Heliotropium megalanthum-Heliotropium longistylum 


0.2362 


0.218 ± 0.015 


4 (4-4) 


1.52 


Heliotropium megalanthum-Heliotropium sinuatum 


0.2498 


0.221 ± 0.010 


9 (9-1 1 ) 


3.42 


Heliotropium longistylum-Heliotropium sinuatum 


0.2288 


0.196 ± 0.012 


7 (5-7) 


2.66 
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H. filifolium/H. floridum 



0.0 0.2 0.4 0.6 0.8 1.0 

H. filifolium/H. megalanthum 



0.0 0.2 0.4 0.6 0.8 1.0 

H. floridum/H. longistylum 



0.0 0.2 0.4 0.6 0.8 

H. floridum/H. sinuatum 



0.0 0.2 0.4 0.6 0.8 1.0 

H. longistylum/H. sinuatum 



0.4 
Fst 



H. filifolium/H. longistylum 



0.2 0.4 0.6 0.8 1.0 

H. filifolium/H. sinuatum 



0.0 0.2 0.4 0.6 0.8 

H. floridum/H. megalanthum 



0.0 0.2 0.4 0.6 

H. longistylum/H. megalanthum 



0.0 0.2 0.4 0.6 0.8 1.0 

H. megalanthum/H. sinuatum 



0.4 0.6 
Fst 



Figure 2. Histograms of F S t frequency 
= distribution for each species pair of sympatric 
1 .0 species of Heliotropium sect. Cochranea in the 
area of Totoral, Chile. 



be under positive selection. This is a signal that allopatric 
speciation may not have taken place among these closely 
related species: Under such a scenario, it is expected that 
a large portion of the genome diverges by a combination 
of selection and genetic drift and that the number of 
divergent portions of the genomes increase with time 
since speciation (Wu 2001; Nosil et al. 2009), which is 
not the case here. Allopatric speciation should also be 
characterized by homogeneous divergence across the gen- 
ome (Wu and Ting 2004). Our data tend to support a 
model of nonallopatric speciation in which divergence 
takes place in the presence of gene flow (Nosil et al. 2009; 
Strasburg et al. 2012). As absence of gene flow appears to 



be restricted to only a few loci, speciation in Heliotropium 
sect. Cochranea is probably an incipient and ongoing pro- 
cess (Wu 2001). Whether speciation with gene flow did 
occur in Heliotropium sect. Cochranea remains to be seen 
when more detailed phylogenetic and population data 
become available, but the data presented in this study are 
at least in line with such a scenario. However, if the dura- 
tion of the periods of allopatry were relatively short and 
populations came relatively early during the speciation 
process into secondary contact, the signature of the DOA 
may be indistinguishable from divergence with gene flow 
(Strasburg et al. 2012), in part because gene flow will 
tend to homogenize polymorphisms that may have 
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accumulated during periods of allopatry in those loci that 
are not subject to selection (Via and West 2008). Luebert 
and Wen (2008) postulated a model of diversification of 
Heliotropium sect. Cochranea in which divergence alter- 
nates with recombination in periods of contraction and 
expansion of population ranges (Stebbins 1952) associated 
with fluctuations of aridity in the Atacama Desert. Drier 
periods would have promoted contraction of population 
ranges and therefore made allopatry more likely. Con- 
versely, more humid periods would have favoured popu- 
lation expansion, thus increasing the probability of 
secondary contact. If pulses of arid and humid periods in 
the Atacama Desert are short and produce the effect 
described by Stebbins (1952), then the predictions of Lue- 
bert and Wen (2008) will be compatible with the data 
presented in this study. 

Climatic and paleoecological studies in the Atacama 
Desert suggest that arid and humid periods have alter- 
nated at different timescales. This occurs as a consequence 
of the El Nino phenomenon and the Interdecadal Pacific 
Oscillation (annual to decadal timescales; Dillon and 
Rundel 1990; Schulz et al. 2012), as well as geographically 
larger scale fluctuations associated with sea surface 
temperatures (millennial timescales; Latorre et al. 2006; 
Placzek et al. 2009; Gayo et al. 2012) and even larger 
timescales related to the glacial periods (Haselton et al. 
2002) or other geological or climatic events such as the 
Andean uplift or the effect of the Humboldt Current 
(Hartley 2003). While most evolutionary studies on Ata- 
cama Desert plants have focused on large-scale fluctua- 
tions and allopatry was therefore the preferred mechanism 
to explain diversification (Gengler-Nowak 2002; Luebert 
and Wen 2008; Luebert et al. 2009; Viruel et al. 2012), 
this study suggests that short-time fluctuations and peri- 
ods of sympatry may also play a role in the evolutionary 
diversification of plants in the Atacama Desert. 

Premating isolation mechanisms might be at play if flo- 
ral morphology is a good indicator of pollination prefer- 
ences in Heliotropium. Floral morphology in Heliotropium 
(Hilger 1987, 1989) is characterized by a conical stigmatic 
head with a basal receptive tissue. A nectar disk is located 
at base of the ovary. The gynoecium is surrounded by five 
stamens in the tube of a infundibuliform corolla, leaving 
limited space to a pollinator to reach the nectar. This 
floral architecture suggests insect pollination, and insect 
pollination has been shown to occur in Heliotropium 
(e.g., Grant 1949; Weiss 1991; McMullen 2007). In Han- 
cornia speciosa Gomes (Apocynaceae), a species with simi- 
lar flower morphology to Heliotropium Darrault and 
Schlindwein (2005) showed that only insects with match- 
ing proboscis lengths are able to pollinate. Corolla tube 
length, relative position of the statements in the tube, and 
ratio between length of the style and length of the stigma 



vary across sympatric species of Heliotropium sect. Coc- 
hranea, but are relatively stable within populations. This 
can be an indication that different species are adapted to 
different pollinators (Luebert 2013). However, whether 
floral morphology has an actual effect on pollination in 
Heliotropium is an issue that needs to be assessed in the 
future. 
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